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ABSTRACT 


A theoretical study of the Replenishment at Sea 
Maneuver, on the basis of mathematical modeling and digital 
Simulation techniques 1S carried out. 

Concepts of Linear Multivariable Systems Theory are 
introduced to implement actual modeling techniques. A 
Mathematical model of a опе plant -two ship system with 
multiple inputs and outputs is developed. 

Risk of collision during the Approach Phase is the 
primary consideration. A thorough analysis of the ships 
natural behaviour under differing initial conditions is 
accompiished, thus establishing points of collision. 

An Automatic Station Keeping Control System, which 
guarantees safety and precision during the maneuver is 


designed and succesfully tested. 
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Replenishment at sea, has become a common and necessary 
practice in almost all the navies of the world. 

Underway replenishment as it is understood today, 
inplies the operation of at least two ships steaming at 
close proximity. Experience and studies over the years have 
shown that when doing so the ships must withstand 
hydrodynamic Forces and Moments that tend to alter their 
dynamics creating situations that involve the risk of 
КОШ слоп. Knowing this fact xt is of pure logic to think of 
some automatic device that could prevent and conpensate this 
hydrodynamic phenomena. Unfortunately, the data availabie 
today do not assure a complete understanding of it. This 
limitation does not mean that such automatic system must be 
seen as impossible; the mathematical model of a skip 
developed by Abkowitz [1], gives the chance to include these 
forces and moments as enviromentai forcing functions. 

Engineer's studies have ranged from the simple case 
where the dynamic of one ship without ezternal forcing 
functions is involved to more complex ones where the 
dynamics of two ships are coupled by the interactive 
effects. 

The most recent one {8 ], which ves used as starting 
το ος this work, treated two Ships Steaming at close 
proximity as one system of multivariable inputs and 
nultivariable outputs, and it deveioped a reliable automatic 
control system that allowed both ships to keep their 
relative stations during the replenishment maneuver. 

This work, investigates the posibility of iwproving the 
above mentioned control system in order to increase its 
capabilities to the general and more desirable situation 
where the automatic control system starts operating at the 


Approach Phase. 
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ДЕ ПАТИПИВЕНАТТСАФ MODEL OF A SURFACE SHIP 


The derivation of a mathematical model representing the 
steering and maneuvering cf a surface ship in open water is 
given by Abkowitz [1]. In its general form, the model 
accounts for non-linear, as well as linear effects. 

Bodies moving ina fluid medium are free to move in six 
degrees of freedom. In order to define the eguations of 
motion, a right hand rectangular coordinate system is 
established, the origin of which is chosen to be in the body 
itself, as shown in Figure II-1. The origin and the axis are 
fixed with respect to the body but movable with respect to 
another system of coordinates axis fixed in space; it is 
assumed that at time t=0 (initial time of the problen) both 


systens coincide. 
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rigid body can be written as two equations -on rorce 


= 


e 
equation and the second a moment equation. The equations 


ае; 


d(monentun) 
(external forces) Fe=— 
t 


d(angular momentur) 
(external moments) Ке=— 


5 


(II-1) 


ПЕСИНЕ БИ τεσπα. ге ship's six degrees of 


freedom have been found [1] to be: 
X=m[ U-RV+QE-X_ ο ον ος τε (PR+0) ] 
= TDYI στι ν ° DUN V προ... p m 
Yan r EI Y (perR2)+Z (PO=P) ] 


Z=m{ W-QU4#PV+¥, (PR~O)+¥, (P+OR)-Z, (Q2+P2) ] 





L=PI, +(I,-I,)OR+m{¥, (W-QU+PV)-Z, (V-PH+RU) ] 
M-QI, *(I, -I, ) PR*n[ Z, (U-RV+QW)-X_ (W-QU+PY) ] 


N=RI, +(1, -I, )PQ*n[X, (V- PW*RU) -Y. (U-RV«QU) ] 


(11 2) 

(where the notation U=dU/dt is used) 
Satisfying the following equations: 

τν > >> | 

Fe=iX+jY+kZ 

— — — — 

Me-iLtjü*kN 
and where the symbols used stand for: 
m Mass of the ship. 
2,142 Components of force in the x, y, z directions. 
Ie M LN Components of moment about the x, Y, Z axis, 
UU, Gomnuonents olsvenoeity іп the x, Y, Z; directions. 
Χς ες «ἲς bastances пеп the center of gravity to the 

ορ σσ the X; y; z directions, 

πο R Components of the angular velocity about the x, 


тахте 
ШІ T Momemts of inertia about the x, ý, Z axis. 


Eta tions II-2 describe the reaction of the rigid body to 
applied forces aS a function of the geometric and physical 
eharseterıstics of the body itself. They do not include any 
of the applied external forces such as propeller thrust, 
rudder forces, moments and forces due to the fins (if any), 
Deectron forces of the fluid (hydrodynamic forces}, and 


Wades and wind forces. 
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Normally, when dealing with steering and 
maneuvering of surface ships in open water, the primary 
motions are considered to take place in the horizontal 
plane, and vertical motions are neglected. As this is also 
the case of interest, where only horizontal motions are of 
concern, the equation in Z, the vertical force equation can 
be dropped. Under the assumption of calm waters, roll, pitch 


and heave are all negligibles, i.e., 

p-P-Q-Q-W-W-0 
Hence equations II-2 reduce to 

X-n[ U-RV-X, ү Е] 

Yy-n[ V*UR*X В-Ү Е) 

N-RI, Fay Z. ъа 0) ^ (UB) i 

(11-3) 

and assuming the coordinate axis origin placed at the center 
ox gravity, Хәй 20, equations II-3 become 

X-n[ Ü-RV 1, surge 

Y=n[ V+UR J, Sway 

N-RI, , yaw 


(II-4) 


The Left hand sides of equations II-4 represent the 
forces and moments along and about the coordinate axes, and 
the right hand sides show the corresponding Gynanic 


reaction. 





O Linearization through Taylor's series expansion 

The forces and moments on the left hand side of 
equations ІІ-1 through II-4 can be expressed as properties 
of the body, properties of the fluid and motion. 

Since steering and maneuvering are of interest, 
forces and moments are also considered as function of rudder 
(control surface) deflections and the change in r.p.m of the 
propeller shaft, furthermore, for a surface ship moving on 
the horizontal plane no forces or moments are due to 


Orientation changes. Hence: 
(X,Y,N)^f (u,r,v,ü,r,V,O6 ,0,5n,etc) 
(II-5) 


In general it is possible to linearize a function 


f(x) by the use of Taylor's series expansion, thus 


E ¿Dx Dix) Ox" Six) 
esta oa tT ES 


where Ax-x-X,. For small values of x the second order terns 
can be neglected and thus considering oniy the following 


expression for f(x): 


f (x) =f (x, ) sx YO 


(II-6) 


The same principle can be applied for small 
Berturbatıons in equations iI-5, which is a function of many 
variables. Since the Taylor expansion is written for a 
particular point, this point is chosen to be an equilibriun 
МЕСОТО: Απ equilibrium position is that of straight ahead 
motion, at constant speed with rudder amidship. The 


hydrodynamic forces and moments have been found to be: 


° Ф & 
= u+ + ^4^X. utX.vtX- rX £ +Y L 
X Х,У X r X. Xe ر‎ 5 „en 


where 4u=u-u, апа Ап-п-п апа ЕНЕ сппесгір 0 іс гетеггей to 


the values of the variables at the initial equilibriun 





Condition and where all the partial derivatives are 
evaluated. For Y and N similar expressions hold. Equating 
the linearized expression for X, Ү апа N with equations II-4 
results in the linearized equations of motion for steering 


and maneuvering: 
X, buteX,v*X reX, U+X. V+X. r+X,Š +X An=mu 
Y. Au+Y v+Y r+Y, Ü+Y; V+Y, F+Y,S +4 An=m (v+ru) 


° ° ° ° 

N,Au*N v*N r*N, U*N; V*N; r*N,$ «N An-I,r 
(II-7) 
,. 5. τ. ην 


vanish for any symmetrical port and starboard shape of ship 


The derivatives X,, X., X 


r u 


(Symmetry about the xz-plane). This has the effect of 
decoupling surge from sway and yaw. 
Furthermore, considering negligible the effect of 


Sn in Y and N equations, equations II-7 become: 
(X; -n)ü*X Au =-X An 
(Y. —m) V+£#Y v+ (Y. —mu) r+VX, r=-Y; $ 


(N,-I,)r*NH r*N, V«N, v =-1,8 


(11-8) 


3. Nondimensionalization 

For computer Simulation purposes, equations II-89 
аге used with the nondimensional coefficients of a ҚКагіпег 
ship, whose characteristics are those of Table JI-2. The 
nondimensional nomenclature is shovn in Table II-1, and the 
nondj mensional] coefficients and conversion factors are shown 
NEU ο πττο ο το it is of interest to point owt 
that the Sca computer tine frane is in 
nondinensionaiized fora. 

In order to sinplify the notation, no special 
symbols are used for the nondimensional quantities, being 


understood that only these quantities are of concern. 


ШЕ 





Symbol 


Ge» 


x? 


Derivative 


Eespect to 


Derivative 


respect to 


Derivative 


respect to 


Derivative 


respect to 


=; 


De ative 
"m 


ла 
Cw 


er 


ес 


(D 


эре 


Derivative 


BDospect to 


Derivative 
кереес το 
Derivative 
espect to 


Derivative 


respect to 


Derivative 
respect to 
Derivative 
PES reet to 


Derivative 


respect to 


TABLE II-i 
NONDIMENSIONAL NOZENCLATURE 


Definition 


of longitudinal force component with 


• ° ° ° 
longitudinal acceleration component u. 


of longitudinal force component with 


longitudinal velocity component u. 


of lateral force component with 


transverse velocity component v. 


of lateral force component with 


<. 


transverse acceleration component 


of late 


+ 


orce component with 


encula cave ls 


wy. A 
yo" CA LÀ 


ὦ 


itv component r. 


of lateral force component with 


yaw angular acceleration component r. 


of lateral force component with 


rudder angle componemt Š . 


of yawing moment component with 


transverse velocity component v. 


of yawing moment component with 


e 
transverse acceleration conponent v. 


of yawing momeLzt component with 
уан angular velocity eomponont r. 
of yawing moment component with 


yaw angular acceleration component 2. 


of yawing nonezt component with 


rudder angle component €. 





πομπεῖα‏ س 


x 


ү? 


m t 


Derivative of longitudinal force conponent with 


respect to change in propeller r.p.a. 


Derivative of lateral force conponent with 


respect to change in propeller r.p.m. 


Derivative of yawing moment component with 


respect to change in propeller г.р.р. 
Yawing angular velocity component. 
Yawing angular acceleration conponent. 


Velocity of origin of body axes relative to 


finid: 


Transverse velocity component of origin of ship 


ахсыс тололу to fluid. 


Transverse acceleration component of ship 


relacije το τα ια. 


Hydrodynamic longitudinal force (positive 


direction forward). 


Hydrodynamic lateral force (positive direction to 


starboard). 
Hass of the ship. 


Velocity of origin Of Ship's aes along the 


X-axis. 


Acceleration of origin of ship's axes along the 


X-axis. 


οπλο Inertia от the Ship «ata respect to ike 


Z-axis. 


Time frame. 





PADLE [1-2 


MARINER CHARACTERISTICS AND MAXIMUM VALUES 
ΝΕ: PDARAMBEIERS 


οποια, Ft 21283 
Bean, ft 10.0 
prart, ft 29.275 
Displacement, tons 16800.0 
Block coefficient, C, 0.6 
Yaw angle (RAD) 026 
Yaw velocity (RAD/SEC) ; 0.0349 
Suge velocity (FI/SEC) 5.068 
Propeller speed (RPM) 30.0 
Rudder deflection (RAD) 05.35 
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ος 11-3 


NOIDTMENSTONAL HYDRODYNAMIC COEFFICIENTS 
NUMERICAL VALUES AND CONVERSION FACTORS 


Nondimensional 


Coefficients 


(Y, -mx ) 
r G 


een‘) 
u 


ШЕ") 


V 


em") 


Note: f= Sea water density 


Nondimensionalizing 


Factors 


0 57pm 
0.5*pL3 
0. 5F*PL*u 
ο. S 
0. 5*pPL*u? 
Vea eo 
0.5*PL3u 
0,520 L2u 
Va 
E AAN 
Do SUD s 
ΕΙ: 
O: SUL? U 
ο ES 
0. 5*Qp*u 
Qus mS 


LAU 


ИВО S Lug Et) 


Nondimensional 


Values*E+05 








4. Computer Simulation 


If the motion of the ship is to be considered under 
external perturbations and with acting controls, no further 
Simplifications in equations II-8 are possible. 

Thrust, rudder and fin forces and moments are 
considered control elements, all other forces and moments 
are not normally controllable inputs, but they must be 
included in cases where the ship has to be controlled in 
their presence. To this category belong the interactive 
forces and moments generated in the case of snips in close 
underway replenishment stations. 

Taking the Laplace Transform of equations II-8 and 
considering the ship steaming in equilibrium conditions 


Me. , Steady forward speed u,=1) 
Ы ш) а J+E(S)[ SY ti-m]--Y1; f(s) 


V(S)[SN; +N, j*r(s)[S(N:.-1;)*H, j]-7N3 dis) 


u(s)[s(X;-n) *X, ] --X nís) 
(177.9) 
S f 
Since r(s)-si(s), equations 11-9 ресопе: 
V (S) [s? (n- Y ) -Y, s JY (s) [-s?Y; *s (a-Y, ) J=¥s (s) 
MUS) [SEN SN, 7+ (5) Е Ен | =} 5) 
ДЕСЕ) (55(4-Х4)-5Х | =X Ὦ (5) 
η 
Metting: 
ΠΥ: 
рі 1=-ү 


v 
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ι1-ῃ 

a^? 1ш-Ү, 
b?1-m-Y, 
сеіс0 
alc=-N, 
bi2--N, 
2 =0 
ET _ N; 
р22--М, 
ce2=0 
ΙΧ. 
b33--X, 


c33=( 


Equations II-10 can be written as: 


v(s) [ al 1524p1 1с+сС11 1- (5) [a21s2+bz2 τσ εἰ 151 ὁ (5) 


$ 


v(s)[al2s2+b12s+c12 ]+Y(s)[a22s?+b22s+c22 Jen ὁ (s) 


ἫΝ 
c | =X, R (5) 
(II-11) 
setting: 
v (5) =A (5) vr 
Y (s) =B (s) Y=B 
u (s) =C (s) u=C 


IFi-Y,Ó(s)-KA1*D1 
IF2-BpÓ(s)-KB1*D1 
IF3=X_ n(s)=KC1*DN1 

Pfauarions 11-11 become: 
al1A+b11)+c11A+a21B+p21B+c21B=1P1 


al2A+bt2Atc! 2A+a22B+b22B4+c22 B=1F2 
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a33C+b3 3C+c33C =IF3 














(II-12) 
or 
3113t*a21B-F1 
Taaa 
a33C -13 
(21215) 
where 
T1=—-p11)A-c11)-p21p-c215+IF1 
I2--b!?A-c!?À-pb2??B-c??B-*4IF2 
I32-b33C-c33C4  IF3 
Solution of equations II-13 yields: 
І1 гез o] ali 11 j ali a21 ΤΊ 
IZ «7^ СО) 212 I2 0 [3112 a22 T2| 
| 
I3 0 a33 0 I3 аз з (9 0 a33| 
л A 775277 
(II-14) 
with 


A = ai2 


22 (0 =а?3 (а11а22-а12аг1) 
3 


and replacing the relations between A, B, C and the original 


Variables v T. u, 
veio v, s [at 
Tere Sipe ate: + [at jet 
Ε O О 
e SO 
u=é=uo+ [Cat 


(11-15) 





The transformation from ship to space coordinate 
system is defined by the following relations, obtained from 


Figure II-1: 
Y=u sinY «v ee 


Х=й боз -V sin Ü 


(mo) 
from which 
НЕСІ 
x=x + iat 
TT 


Equations τι; through II-17 were translated into 
WoL 360 Conputer progran: With a constant rudder 
Reflection 5-р1-0.1, the results are shown in Figure 11-2 
(yaw dalyic Versus time; and Figure II-3 (suav versus surqe)l. 
ime Characteristic turning radius of the ship. It must he 
Considered that the described trajectory is valid only for 
small rudder angles, since linear theory does not give the 


Speed reduction in the turn [i1]. 


Sr Controlar rana stability Ol. спе linear model 


Owe oar aner amt qme o مھ‎ Ome Dar aar TETO سیت سید س ست‎ mens w — ane < < cu us eer s F 


These two concepts form a basic hackground to the 
subject of the qualitative characteristics in the handling 
Of ships, since each one takes into account the intrinsic 
ШИӘрегізес ой а Ship and the transfer of initial state of 


ο που tc another state in a finite time {10}. 


ПРЕ ОЙ ЕСО атату CII 


Considering the systen 
UA (t), vlt] for t2t,, and x(t,)=x, 


(11-18) 
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If there is a finite time t 2t and a control 
ut), telt,,t, |. which transfers the state X to the origin 


at time t, + the state x is said to be controllable at tine 


[i - 
O 

If all values of x are controllable for all to 7 
the system is completely controllable, Or Simply: 
controllable. 


Kalman í 6 ] has shown that a linear, time 
invariant system is controllable if and only if the n*ın 


matrix 


(II-19) 


has rank n. If there is only one control input (п-1), а 
necessary and sufficient condition for controllability is 
that n*sun matrix E be nonsingular. 

Taking equations 11-8 without including the surge 


equation (because it is decoupled from yaw and s 


зг, 


„4 
MN 


equations) and dropping out the terms containing the 
hydrodynamic derivatives Y, and WN; in the yaw and sway 
equations [29 ,Pag. 472] only to apply the Kalman Criterion, 


the state equations are: 


e 


-bii/aii -p21/g11 ү y 


toad 
` ç 


-р12/аг2 -рг2/аг2 Js N 
(11-20) 
after performing the arithmetic operations 


zug OS 


ES 
' 
1 


270012 6085 


which is rank 2, indicating that complete controllability is 


assured. 
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D Stability 

The stability test determines whether or not the 
ship returns to an established equilibrium condition 
(straight ahead motion at constant speed), after removing 
the small disturbance which caused its departure from that 
equilipriun. 

A dynamically unstable ship cannot maintain 
straight line motion when the rudder is amidsnip. 

The behaviour of the ship can be analized by 
considering either some introduced disturbance and zero 
control input (4-0) ог the control acting as a disturbance. 

For the first case, and neglecting the surge 


equation, equations II-8 reduce to 
у (5) (а115+р11) +г (5) {а215+р21) =0 
v (s) (аї 25+Ь?}2) +г (5) (а? 25+р22) =0 
(11-21) 
yielding the characteristic eguation 


alis“+p11 altis+h?1 


al2s+b12 ассегрее 
OL 


(a!*a?2?2-ai!2a?1)s24(a!1b224a22b11-a12521-a21p12) s 


+ (b11p22--p12p21)=0 
replacing values and rearranging 


s2+0.685s+1.016=0 


ο ο os S belong to the left half s-plane; the ship has 


fixed control stability with characteristics [1] 





cO. =1.008 
£ =0.34 
Figure II-4 shows the Root Locus corresponding to 
lion 11-22, which is obtained with Fortran computer 


program 11. 


6. The transfer functions 
Defining 
Kil=qglis2¢bliis 
K2l=a21s2+rb2ls 
ІР 3І2гегірігсе 
К22=а22с2+}2 25 


КзЗ=аззсё+рз зс 
equations II-11 can he written as 


v (s) K!1 (s) K21-Y d (s) 














u (s) K33 =X n(s) 
(IE) 
J 
СОИ о па Lor у (5), V ic) cuts) 
νυ ται ο K11 Yrd(s) O K11 K21 Ya ó (e) | 
wed(s) K21 0 к1г wgd(s) 0 Kiz K22 NpO(s)| 
24 TS} O += τη πη "m 0 0 A mss 
nl o o — u р а оссо тт -- 
(11-28) 
еһсге 


eit қ23 : 


= Kee ο 


9 О КЗЗ 


------ 





and replacing the K's 
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Ass (a33+b33)[ (s? (а11а22-а1 2а21) ts (ai!b?2*a??p! 1 


-al2p12-321p21)* (b11p22-p12pD21) ] 
Evaluating the solutions defined by equations I1-24 


u(s) Ku 


п (5) stt 








у (5) Ку(5%2у) 


S(s) s2+ps+q 


V (s) η 
d(s) s(s2+ps+q) 











(11725) 
where 
атая 
5 —  —rrp, — TA 
а21а22-а12а21 
Nga! лү, а12 
r= <. ә وہ 4 ص‎ 
αν Σα στα 2 да? + 
X 
Ku=— 
a 33 
Ys b22-N, b?1 
ае 2222. 
J; a=; a3 
Ns bi DES 
Baer aa == s= 
Ny ati-yyate2 
ὮὉ 3 3 
> 
a 33 
eae oe ey hte ag hehe lan 2ij, ie 
P mm m m —— nn 
allqgee-qgie2gel 
Peet ye a=) teh < i 
alla22-a412421 
(ISO GE 


The transfer functions defined by equations 11-25 
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and II-26, together with the coordinate transformation given 
by equations II-16 and ІІ-17 lead to the block diagraa 
representation of the ship, Figure II-5. 

The numerical values for the transfer functions are 
those of table II-4. 


TABLE II-U 


NUMERICAL VALUES FOR THE TRANSFER FUNCTIONS 


Kv- 0.21447 
Kr=-1.91507 
Ku= 0.00588 
zu 23270923 
zE= 21.227369 


p= 0.68467 
g= 1.01659 
t= 0.18117 


БО SHIPS UNDER FORCE FIELD EFFECT 


The Simple case where the equations of notion of one ship 
were derived (ІІ-9) did not take in consideration any 
External field forcing function. 

As it has been previously stated, when two ships steam at 
close proximity, the existence of this Venturi type effect 
ЕСІ ПО longer be ignored [1]. 

Gaivano [a] using available information and digital 
computer work synthesized the Force and Moment field in a 
set of curves, which can be considered as valuable 
information in the designing of an automatic station keeping 
control system. These curves are shown in Figures I1I-6 ana 
II-7 respectively. 

ita o dare culty in handling these curves is the fact 


ка eye annot be written in a simple mathematical forn, 
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but the availability of digital computers minimizes this 
difficulty to the point where they can be stored as Discrete 
preknown information (Look-up Table), this is shown in 
DSL/360 progran III, Subroutine Slopes. 

А11 discussion concerning interaction hydrodynamic 
phenomena, so far, has been for the Mariner which is assumed 
to be ship A of the two ship system,in Figure II-8. 

But hydrodynamic derivatives for both ships, at any 
relative position of interest are needed to write the 
equation of motion for the multivariable system. If the two 
Ships are identical the force or moment felt by ship A at a 
given Dx and Dy are the negatives of the force or moment 
felt by ship B at the same Dy at Dx=-Dx [4 }. 


Defining 


Y=rorce in ship A due to ship B 
Y=Force in ship B due to ship A 


N-Moment in ship A due to ship B 


NW=Homent an ship В айе 1 hip A 


LO Ship 
then 
Y =-Y 
D z Ot D x - — 5ο 
р Оу : В 
апа 
К =—N 
Ох = Dı = — < 
Dy = fo Dy = b 


(11-27) 


where & and © represent numerical values in the range of 


| 


interest of the variables Dx and Dy. 
The equations of motion of the two shirs, coupled now by 
the interactive effects can be written as follows: 


: {17 X ) 
v (s) [al 1 5Ş2 + þ3 Les ti 1%11(5) [a?*is?tb? 1er 91 сүү d (s) RUM LEN, 


1 
=; 
Y u : - 
v (s) (а1і252%0125%с12 CS) fae=s-+pe- 264022 ]=1;9,05) *N(s) 
$ 


у2(5) [ al1is2+hb11s+c11 ts) Гагісетрг 1с̧+с2 +} ]* Y; f(s) +Y,(S) 
3 





MS [a12s2+b12s+c12 ΤΙ. (а2252%р226%с22 ]- 55 dts) tN (S) 


$ 


u,(s) [a33s2+b33s+c33] ins) 
1-(5) Газ352%0335%с33) XS) 


$ 


(II-28) 


in which the coupling effect is implicitly defined py 


equations II-27. 


1. Computer simulation 

So far, the equations of motion of two ships 
coupled by the interactive effects (II-28) are in a neat 
form to be simulated in the digital computer, but it is 
desirable to impose an additional constraint to reflect a 
real limitation. There exists a certain time lag between the 
mutant the rudder order is given until the action is 
completed. This time Lag 5 taken as σα (non 
dimensionalized time, actual value is 2.08 seconds}. Thus 
the rudder order fis modified to take into account this 


time lag. 


dis) NEU. de 


oc 


TM 
(11-29) 
where 
de is the desired rudder angle 
б is the actual rudder angle 
tr is the tine lag 
a. The egustioP of motion of ship # 1 (leading ship) 


From esuations 11-28 


.Τ-- στον SG) αυ [ς δ (5) HESS, 


і 
$ 


7 => - 7 2 5 > = c З E? 
v (s) fai 2524р1254сі12 Ls) (а2%252%р2 eet Cad Ј= %,, d (s) 4 ШІ) 


% 
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м р аззаг+рззс+сзз ] =X, n (S) 


(11-30) 
Considering the steering control as T Ë is 
described by equation II-29 
_ 10 0,,(8) 
5%10 
(II-31) 
and setting 
Υεόμι (5) 
ο *T(5) 
0.1s+ 
IF21 (5) _ τ +S) 
K 0. 1s+1 E 
(II-32) 


The same procedure followed for the derivation of 


il-12 through can be applied yielding as before 


e 5 


te е 
ΤΟ a2 D218 +C2!B, =IP11 


ος % ee e 
EE Reno CLEA 134558 але ccc p-rEZ2! 


a33C, *b33C, «c3 3C, -IF31 
πι a) 
or with 
1112-bi1i CIR -b21B,-0228 «ІРІ 
I21=-bi2A -cı2A -b22B -c22B +IP21 
131=-b33C -c33c, +1F31 
(а 


equations II-33 can be written as 
te 96 
ah жа<%8 -I11 


al2A +a22B =121 
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азс. =131 


(11-35) 
Solving for A, B, and d 
‘e а22111-а21121 .. а11121-а12111 .. 131 
(II-36) 
and the original variables are recovered as before, 
Y, =A, =¥, + A, dt 
Y, -5, +B, кф +] B. ‚|>, dt jdt 
Ὁ τὸ τῇ Ἡ Ὁ, ας 
(11-37) 
so that in the space coordinate system 
Y =u, Sin, tv, Cos Y, 
X, zu. соз? =V Sin is 
Y =y tli at 
x, 2X, e| X, at 
(II-38) 


DEED C C OD οἱ ος motion ot ship $2 (trailing ship! 

Since, from equations 11-283, at is seen that 
equations for ship #2 are the same as for ship #1 differing 
ШЕШУ іп the subscripts, this derivation is omitted. 

HUMO Io Oro Cs ap Iv ana. Ship’ $2 were 
Mas tatea iato DSL/360 digital computer program IIZI. 

A piecewise linear approximation of forces and 
moments is given py the table look-up and the interpolation 
Subroutine Slopes. A warning message is printed whenever the 
distance between the ships become less than 25 feet. For 


separations greater than 250 feet the ships were considered 
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to pe outside of the range of interest and the forces and 


moments assumed to be null. 


AS previously mentioned, the existence of risk of 
collision in the replenishment maneuver is qualitatively 
p known [3] and it fixes the basis for any attempt of 
designing an automatic station keeping control system. 
However a well founded quantitative analysis of this 
particular subject is almost impossible to find in the 
current literature. 

The open ioop test presents a quantitative analysis 
of the two ship system behaviour, and by choice of different 
Mal.” conditions it detects the Collision Points. Of 
course, the simulated situation is ideal and thus only 
approximates actual physical conditions. The analysis is, 


however, an excellent approximation to reality. 


Шы ШИ шу ао 15 ideal іп the sense tha 1 
intervention 1s excluded and only the inherent 
@mmeacteristics of the ships іп the hydrodynamic field are 


evaluated. 
Table ІІ-5, indicates the tests performed to the 
multivariable open loop system using 50 рата 


Computer Program Ill. 


ТАБЕЕВ στο 


TeSt i U1(0) U2 (0) DY (0) DY (6) 
1 m ТЕ 622 - 1. 
2 1. ie Шс? 0. 
quo s Dos O. 
Ц 1; ive 2 =a 
5 12 1.2 e O. 








Note that the actual distance 1S obtained by 
multiplying the nondimensional length by the length of the 
ship (i.e., 0.2 <=> 0.2*528 feet). 

Figures  II-9 through II-20 show the behaviour of 
the two ships for the conditions pointed out in Table II-5. 

Three important conclusions are obtained from the 
analysis of the Figures. 

1) The ideal case, where the equilibriun conditions 
are set at abeam position at t=0 and released at TON M 
= the combined effect of Force and Moment fields tend 
to pull the ships apart. It is apparent that no collision 
Mier is involved (Test 2, 3, Figures II-12, II-13, II-14, 
II-15). 

ii) The more realistic case uhere the equilibrium 
conditions are set at the Approach Phase and then released 
ο), the combined effect of Force and Moment fields tend 
to attract both ships and collision is unavoidable. (it must 
be considered that after t=0 no human action is involved and 
οσους observed in Test i and 4, Figures II-9, II-10, 11-12 
ИЕ иге II-16, 11-17, {1-18 correspond to the natural 
behaviour of the ships). 

NUES το ο tigures-II-19 anc TT=20, Shous the 
Departure Phase, where the speed of ship #2 is allowed to he 


uucrGased and no risk of collision is observed. 








1119 THE AUTOMATIC CONTROL SYSTEM 

This section starts with the definition of a plant 
With multiple inputs and multiple outputs and proceeds with 
mathematical derivations until a model for the automatic 


control system is obtained. 


DE Ur MULTIPLE INPUTS, MULTIPLE OUTPUTS PLANT 


The initial equilibrium condition for two ships steaning 
at close proximity indicates that an undetermined d. forcing 
function is required at со to compensate the effects of 
F (AX? AY, ) aná (AX, + AY) which are considered constant. 
This definition is valid only at t=0 since there is not any 
logic that prevents ó from changing ὁ, το another value that 
compensateS increments of F and lh due to variations of АХ 
and AY. Abkowitz [1] has shown that the force or moment 


- 


caused by change in any var 


[45 


able is expressed as the oroduct 

of the derivative of the force with respect to that variable 

(with all other variables at equilibrium values) and the 
К 

change in tne variable. Using this criteria, at t=0 it is 


required to compensate: 


ЭЕ Р 
P (AX, AY) =o AX *2—AY 


ox Oy 
and 


` 


ы 





“Oh 
M (AX, AY) =—AX+-—ay 
OX 


оу 
(ra 


which represent the "departure" from equilibrium conditions. 


Defining 
“OF (AX, GY) =ki9 AY+k20Ax 


and 
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"MM (AX, AY) =9*0AY+q?0AX 
(III-2) 


where k10, k?0, q10, q?0, represent the rate of change of Е 
and M with respect to AX and AY measured at AX, and AY, «the 
longitudinal and lateral separations from the established 
equilibrium conditions. 

The linear expression for OF and OM must be defined in 
terms of the variables u, у, Ф апа then added to the left 
hand side of equations τ) under the following 
assumptions: 

a. The two ships are identical. 

b. All hydrodynamic coefficients are not affected by the 
intermingling of the water pressure between the ships, 
therefore remaining constant. 

C. The two ships are considered already as being 
Eongside edch other (AX,-0). 

а. The forces aná moments acting on the ships are equal 
ШОШ ООПТ йс and of opposite sign. 

e. The change in forward velocity is considered 
negligibie. 

The increments AX and AY can be expressed as given by 


equations II-17 
Ax-x, -x, - |d, -X pat 
r Y VY =Y, =| (i, zu nue 


Using equations II-16 


" 


Шз» .... А V) E “ Vr) XT 2 RM : 12 
X, ^ Eau. Cosi, ~v, Sin -t, ) (U, Cosi, М, sind, ) 


F 
% жу, сов?) 


= ) [7 ? cM - 1 17 
Y, Eu. = (u, sin +1, Сос) (u, San 


(1ЈІ-3) 
For the small perturbations being considered 


Cos T «4 sin? ~0 





Е D 
f 


cost, 531 Sinf, ^0 


and equations III-3 reduce to 


= — ma 
Y, Y V, Y, 


and finally 


AX=0 E ΔΧ(5)90 
Are) (vy -v, )dt L ЛҮ (5) =—(v, =v) 


The equations II-10 are modified to include the 


interaction effects and become 


а е Ие, Js) Pas DS (Ss) ae (s) 


O ІҒІ(в)(а22524р225жс22 1%У:(6) 4-8; d (s) 
$ 5 
Vo{S) [a1 σοι 1ςσ.οί 7505) [ a? 152-:521sgtc2i Ј+ (s) K=Ys d S] 
$ $ 

^ 
ts) ING OS) 


——© 2 
$ 


1 
J 


a 
N 


£m? 
= 


ate 
«^ 


en A . Er `y 2 
w(s) [a12s2«b1?s«ci? УД (з) гаг2521р2 2: 
5 | 


(111-4) 
where now 

cil=-k and ci2#=-q 

letting 


р%9-а1:524р1 "στο. 
р20=а21524+р2ї5+с21 
p30= merece Pyilicatcie 


p*0=a22s2+h225+c22 
(ІТІ-5) 
equations 111-4% becone 
104457, c 20 f JR 
(5) p1 °+1,(s) p2 °+ va(5) k=Yg 0(5) 
> 


vı(s) p39 4 (s) p+0+v(s) q=ty d,(s) 
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5) K+ (5) p1o« Lis) p20-1, J,(s) 
jS) ч+%(5) рзо+Ї (5) р+0= 9 (5) 


(ILI-6) 


Баага опо —III-6 Show that two ships affected by 
interaction forces and moments can be described as a 
multivariable system where the deflection of the rudders d, 
and d, are the control inputs and the yaw angles Y ana Y the 
outputs of interest. The next following step has as object 
to determine the form of the entries of the open loop 
transfer function matrix G (whose model is shown in Figure 
III-1), so that the system can be analyzed and modified, if 
necessary, to become steady state decoupled -after a 
transient period of time, a variation introduced in the 
rudder angie of one ship vill not alter the yaw angle of the 
other | | 


Eron figure III-1 


η 11 22 ) 
d(s) = | ; E 3 3 s S +, 
/№ 


οἱ У ЕО 
2 Y 
Or 
T T 
Ὧ (s) - dts) 6 (s) 
NS 
[1 
AT 

Solving for Y (s) 

Ç 7. ы 

Od dg σι: 

д2 g<e 

hence 


Er. =0 91 1+2: 
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= азге Догг 











(III-8) 
and evaluating the gains 
Sa р 
di 4-0 1 jo 
E Ы 
Í "η 
1 19=0 1 = 0 
ES IS) nust be of the type 
ho B 
d1 di 
G(s)- 
y, 9; 
d> б 
τι 
1. Trapsfor function matrix type number 
Бастон 11-53 σαντα the general form of the Plant 
EI ο Іп explicit expression for the entries of the G 


Matrix is not needed. The steady 


States by 


knowledge of the number of pure integrators in the 

of G, this is referred as the type number of the Ma eb & 
: ] 

Solving III-6 for Y,(s) 


6 


Tail S 21 Ji 
E cs z Y 


f 


S) ΜΞ ο 5] ]Ἰ 
ele 


i A 














where 


ОН ТЕ e) and 


State decoupling of the 


means of a compensator matrix GS requires only the 


ent ros 


i 5 je 
3 (2 ` r 
recaltimg 


(111-10) 


(III-11) 





A = (p10p< 0-p20p30)-(qp?20-kp+0) 


(EX 1-12) 
where рі0 and p39 are second order polynomials ins 
satisfying 
Lin p#0 
s->0 


mus the special case of having either k or q 
identically zero is avoided. 
solving for Y (s) in III-6 


417 (Ns pi9-Ye р39) (р19р%0-р20р309) f 
Te di) 


= ( 40-. 203 
4 De q- a 9) (kp*0-9p20) $ (a 


αν 


A 


(111-13) 


where A is given by equation TT I=12. 
Replacing p19, pee, p39, p*® as indicated by 


equations III-5 and taking separately σα11 (5) and g??(s) 


us) (К ОО 5-17 (ас25верегвясег) 
qi 5) See a ns pe 


ὁ {5) Ж 


Кот (5 ΙΙ с 1) (ағ21524р2 Larc? 1) (a? 2024 B1 estei?) 


(opcm s vm s 00-2 


Үс (а1252+рі lstci 1) (aleésetblesicie) (a? es?+b??s+rc?®) 


ص 


А 


Aum مو‎ 





К €— P ЧР 

















w Geer sen a erie 





ا — .< .— 


38 








νι ο ο (C J) (асері 25жсісужҮ;ак(ае2сетрг2сжс>2 2) 
а : Т 


 Ү42(а2152%0215%с21) «Қх2 (а2252%р225%с22)-Қақ(а2152%р2154с21) 


Д 


and from equations II-10, II-11 


(III-14) 


ο 22=0 


then one s can be factored out immediately in the numerator. 
The independent tern must be also checked, if it is equal to 
zero, then another Ss can be factored and So on. 

ер ас С Шу =F, and, с12 by -а, the 


wh 


independent term becones: 
Ngk2b22-Npkqb21-Yekqb22tYe b21q24Ye gkb22 
=Y9q2b21-N9k?*2b22+Ngqkb?21=0 


"ch 1Jnaicates thet @ second S can be factored out іп tho 
numerator. 


Expanding the denominator 
A =s2[ (as3+bs2+cs+df- (est+dy | 
(\=s53[ a2s5+2abs + (b2+2ac) s3%2 (adtbc} s2 
mco:-2bd-e)se-2dí(c-e) 1 
where 


a=sal!a22-a2igl2 
b=bita2Z2¢aqlip22~q2ipi2 
Czbilb22-ka?24a21q-b21b12 
d=b2ig-kb2?2 


есда2:-Қаг2 


It can be shown that the independent term never 
goes to zero for the given values Gf the hydrodynamic 


coefficient of the Mariner [9]. 





Hence, the highest factorable power of s is s3, and 
gil ts) 93209) 
1 (5) 


is a type 1 transfer function (i.e., it has one single pole 
at zero). 


Expanding the numerator of g21(s) 

Ма4а21 (s) =s[ Ye k (a?*s+b22) (al2s2+bl2s+cl2) 
—kN£ (a21s+b21) (a12s2+b12s+c12) 

+Nsq (alls2+b11s+c11) (a21s+b21) ] 


and making the sane analysis done before for the independent 


Bern: 
SEDES Yare tin DSlg-Npqkbel=0 


Which indicates that the highest factorable power of s is 
ESHence g?!(s) is a type 1 transfer function. 


е7 


ne Dom or eguation III-1! for Lí(s) xs 


V 


(Ne p1*0-Yç p30) Er DD 
— s 


(Ys g-Ng k) (kp4t0-gp20) 1 
ол 


say з, у | 
с == GI 








^ 
/\ 
= 
ο το 
a= Could be expected by symnetry, το has the same forn as 
ТУ (в). Da conclusion, the type number matrix is in 
general 
t 1 
Tp= 
AY 
Ша 71 
(ТтТ=1р) 
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2. Steady state decoupling 

Once the type number of the plant matrix G has been 
found, interest is concentrated in the determination of the 
cascade compensator matrix. 

An important consideration that must be made is the 
fact tnat the technigue to be used in finding the 
compensator matrix does not assure by itself stability or 
good transient response when the feedback path is closed. 
Figure III-2 shows the closed loop block diagram, where Сс 
is the conpensator transfer function and Gp IS еа руапЕ 


transfer function, and 


dis) Yıs) 
R (s) = . «(85 


ЖЕ Y 


2\5) 
Huang and Thaler [5], had shown that assuming that 
Gp is stabilized through the eOnT 2gULraLLoR.orf Trigu rinna 


then the system is steady state decoupled if and only if 


Li (3. $6CpGc]i, 


lim — ern 
E) S det (Z*GpGc) 


for all i-1,---------,n 
(iII-17) 


If Ge is a diagonal matrix and all inputs are steps 
{kj=1), the conpensator type matrix TC, which gives the 
nunbe of integrators required for steady state decoupling 


" 
ОЕ а 2%2 system is obtained by satisfying tiie conditions: 


(111-18) 





where 
M2Max[ (tpll+tcC11), (tp+2+tc22) ,(tp!i+tp22+tcl1i+tc22), 
(tp12+tp21+tc1l1+tc22) ) 
N12Étc114+tp21 


N21@tc22+tp12 


(711219) 
Eod equation III-16 
ΤΙ 1-51 
then 
N12=tc11+1 
N2i-tc22-41 
M=Max[ (tc11+1) , (tc22+1), (2*+*tcll+tc%#2) | 
Considering that tne plant has TOLL pure 


integrators, no more are required in the compensator matrix 


then 
tctt=tcee=0 
and 
Ni2=N21i= 7 
H= 2 


Satisfying the conditions established by equations III-18 


From which: 


rj 
N 
H 
tu 
Q2 
Q 
И 


t 
° 


[11-70 


Επ πο To and vith 011-20, 6 will have the sane 
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type number matrix as Gp, hence 


Mme 


T= 
AS 


then G can be written as 


P12 (5) p12 (s) 
SPA(S) SPA (S) 
G= 
^h, 
Рё! (5) EOM 
SPA(S) SPA (S) 


Where Pals) and Pi, (s), i,j-1,2 are polynomiais in s such 


- 13] 


that 


r (o 


Ss“ 


Ej 
O 


το τους 
de NA Ed. e - 


ὃ 
> 


aa 


means of equation ii1-18 


— 


Pi!(s)tsPZ (S) P32 (s) +sPA (S) 
SPA (S) «ο SPA (5) 

р21 (5) t sme (S) P?? (s) *sPA(S) 
SPA (s) SPa(S) 


(pilgs) *sPA(S)) (P** (5) *sPA(S 
Ор ту с s Ga nun 
p (sPa (s) ) 





-s 





-—— 


(P21 (s) *Pa (5) ) (P12 (s) *sPA (5) ) 


(SPa (Ss) ) 





= = 





E Pel Da ο) 
E 
M SB (s) 


D cus ИБ is determined by 








р12 (5) +5рд (5) 


Е = 
2+6, SPA (s) 


1 [2+6 ]12 1 [I+S la 


lin КОЛЛ со Кооз 9 
5->0| 577 det LI+G] 991-0 det( I*G] 


and after simplifications 


$255, (s) ME p12 (5) 
lim uu eee ccc eme m) 
S-»0 Pll(s)P??(s)-P?!'(s)P!* (s) 
ο K2 52} 


(III-21) 


Then the system 1S steady state decoupled not only 


NUT Step inputs, but also for ramp (k=2) inputs. 


ОО wir CLOSED LOOP SYSTEM 


Theoretically, the steady state decoupling of the outputs 
of the plant matrix has been guaranteed and the order of the 
system kept low because no extra integrators were added in 
the compensator matrix Gc. However there stiil remains the 
problem of Satisfying an adequate performance, which in turn 
1s the main purpose of the automatic control systen. 

іва Гв] solved this problem by modifying the compensator 


EStrix to the forn: 


δη 0 K1+sKT1 0 


Q 922 (5+222) 0 σου το 
ο — 


and satisfied the desired performance, which was initially 
stated in section A. 
στα exhaustive analysis of the control system, the 


idea of increasing its capabilities based on the conclusions 


ац 





coming next was a feasible and attractive prospect. 

i. Existing decoupling in the Surge equations of the 
multivariable system allowed minor changes in their forward 
Speeds without appreciable effect in their trajectories. 

jy Equations II-27 indicated that no additional 
ШЕТІСПІЗСісс were introduced by changing initial conditions 
the X axis (i.e.,40X*0). 

iii. The gains of the compensator matrix obtained by 
Parameter Optimization Techniques [ 8 ] allowed different 
damping rates in their transient responses. 

lv. No automatic speed control was required, since minor 
changes in speed of the trailing ship vere easily made by 
ОООО т ПО the acceleration tables of the ship (i.e., it is 
assumed that the leading ship forward speed is kept 


constant). 


j. Generalization of the multi: 


(m 
|" 
y 
| 
lo 
2 
[{η 
{ΓΕ 
ТЕ 
th 
Ф 
ζω 
γω 
3 
(b 
FL 


— —— mee = = с зш с шз шыш 


—— w wo х= =ч Den nn 


Only a ¡ey words are enough to describe this 
Objective. Reliable work of the automatic control system, 
capable of keeping the desired relative positions between 
two ships from the beginning of the replenishment maneuver 
(Approach Phase) to the end of it (Departure Phase), 
Em zing the risk of collision. figure III-3 shows the 
general case where both ships nanevuver. 

ο ο ο πο ο GSE interest in the control loop are 
г ᾽ 6 s 
mU. y. у DY, DY. Although DX and DX are also 


Variables that must be controlled, they do not affect in tae 


5 loop, tnezr action is controlled manually hy 
changing the forward speed of ship #2 until DX is reduced to 
Ero. 


The design procedure is an extension of the work 
presented by Lira [в], in which the Ships were initially 
considered steaning alongside at a lateral separation 


greater than that required to perform the replenishment 


τ... 
I 





Maneuver. Tae optimal gains as given in Table III-1 were 
obtained as result of choosing the desired trajectories for 
each ship and defining a cost function to be minimized. 
These values were kept the same for these studies, however, 
the initial conditions, specifically U20, DX, and Y1 and Y2 


were changed in order to include tne Approach Phase. 
2. The controlled pliant equations 


The block diagram of the controlied plant is shown 


in figure III-4, and the following relations are obtained: 


f 1 [(KTY1S+KY1) Y1(S) +(KTls+K1) Š (s) ] 
nz i 





trs+1 
(TIT) 


1 [KY2 (DY-DFIN)*KTY2SDY* (KTZs*K2) Y,(s) 1 





4 (з) = 


42 trs+1 
(III-2U) 
where DFIN is the Desired Final Lateral Separation. 
ποτ ο οἳ ева there exists limitations in the 


maximum values of the rudder angles [3] imposed an 


additional restriction, which is 
44130.359 or (X20 ); for 14120. 


11-23 аһа 11-24 ресопе 


d (== στο * (0.349) 
MMS 
d (s P. * (0.349) 
e trs 1 
(11-26) 
Digital computer program #4, which is progran ¢3 


КОБО ООШ оста тла equations 1123-23 through I11-25 allows 
us to make the simulation of tne multivariable system now 


implemented by the station keeping automatic control systen. 


б 





TABLE- IIL- 


FEEDBACK LOCP GAINS 


Parameter Run #1 Run #2 
U10 de m 
U20 122 2 
DY (0) 0.4 0.36 
DX (0) -1. -1. 
DFIN 0.24 (0% 3 
Κ 1 24913 ο οσο. 
KT1 Dans 3.080 
K2 ΕΠ 2085) 
KT2 2.042 2.141 
KY 1 3.008 3320 
KTY1 2.703 55 
КҮ2 ie t 582 
СТЛ 22818 25600 
f gures ITI S through III-10 Show the computer 
Simulation results of the controlled plant. Both ships 


reached the desired steady trajectories allowing a safe, 


well 
Separation error oxists 
be 


as precise, 


which can easiiy 


naneuver, 
(C202 <=> 


removed by indexing in the hardware. 


although а 


approximately 


as 


residual lateral 


11 teet}; 


Wome that these are not solely I. C. curves but are tests oi 


System behaviour as 


speed 20% greater than that of ship 


Seeeursion of ship 1 is 


caused by the pressure field, 


ШИ al conditions. 


due 


T The 


and in part by 


<. 


i; 


initial 


Ship 2 approaches ship 1 at a forwara 


Par e 


POAN E OEL large Nolen 


23 ” Y ρον 
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IVe A NEW OPTINIZATION CRITERIA 

The preceeding section led to the final design of the 
automatic control system. A basic assumption was made in the 
sense that the feedback gains were kept unchanged from the 
optimals obtained by Lima [8]. 

Computer Simulation demonstrated that the assumption 
made was good enough to allow the two ship system to satisfy 
the desired performance. However tne trajectories can no 
longer be said to be optimal, since completely different 
initial conditions were imposed. To overcome this fact two 
techniques were considered -~the first one, to reformulate 
m optimization study incorporating the new initial 
conditions, and, the second, to use the actual gains and to 
@eraginate a systematic procedure, basically Trial and Error, 
to determine not optimal, but guasi-optimal gains. 

Because it was expected that the range ef variation cf 


the parameters would be reasonably smail, the second 


Fh 


procedure was adopted. O course, the @riel and ЕГЕСІ 
procedure has the disadvantage that many statist 

combinations are ignored, put, on tne other hand experience 
gives a high degree of confidence in a systematic computer 


trial procedure. 


EENEURIAL AND ERROR PROCEDURE FOR QUASI-OPTIRIZATION OF THE 

GAIN PARAHETLRES 

Selecting the time basis trajectories shown by Figures 
ШО and Ili-9 the new objective 55 to smoektha Suay and Yav 
transient responses. Note that the sane grocedure applies 
EN Ugures IiII-5 and III-6. 

From Figure III-4, a Procedure Flow Chart was made which 
includes all the possibilities of discrete changes in the 
K's parameters in the range 0 to 100% by increasing its 


values. 


58 





The possibility of decreasing the K's vaiues was left out 
because it was estimated, pased on linear theory, that it 
would not meet the requirements. 

Figure IV-1 shows the Trial Procedure Flow Chart. 

From computer simulation, many runs were done using 
DSL/360 computer program #5 but after the selection of the 
"best set", graphs of Run #1 through Run #4 were considered 
important enough to be included for purpose of analysis. 
Note that ali of these Simulation runs were made with the 
trailing ship approaching the leading ship at a speed 20% 


greater than that of the leading ship. 


Demet kL AL AND ERROR PROCEDURE RESULTS 


Table IV-1 shows the gain parameters for the four main 
runs of the Trial Procedure Flow Chart. Figures IV-2 through 
IV-9 are the computer output for this stepped procedure. 

ММ) 12:52:12: СІ tho best set of K's was somewhat 
Subjective and based principally on the shape of the 
corresponding computer output CUEVES. A yualitative 
analysis o£ the Sway and Yaw curves, however, indicates that 
there is a definitely improvement in the smoothness of tne 
E UuSrent, Figures 1Y-8 and IV-9 being considered the best 
K πα be seen by comparison with Figures 11I-8, ILI-9. 

Figures IV-8 and IV-S are the most important because they 
assure that the systematic Trial and Error procedure vas ап 
excellent “τ. (Ппиав!-орыкілаі gains for the 
specific initial conditions of the Generalized Desired 


Er Uornance described in section III. 





ΤΕ ШІ) 


FEEDBACK LOOP GAINS FROM TRIAL PROCEDURE 


Parameter Run #1 Run #2 Run #3 Run #4 
KT 1+50% KTY1+50% EY 2m K1+15% 

U10 12 le T: 1. 

U20 1722 {ΠΠ 122 1 

DY (0) 0. 36 0.36 0.36 0.36 

DX (0) -1. -1. -1. =e 

DFIN 0.3 Ca 053 0.3 

K1 3.069 3.069 3.069 2502985 

KT 1 4.62 {Ὁ ὁ 5452 4.62 

κ. 25853 2.553 2:853 2 05: 

KT2 2 141 ӘСІП 2.141 2.141 

KYI 3.320 3% 520 2220 22229 

КТҮ1 2299) ο) 3.8325 5263525 

pr 93 ος” Τ᾽ 721075 

KTY2 2.608 22008 2-606 255 





Ve CONCLUSIONS 


It has been demonstrated that the intrinsic dynamics of 
ships steaming at close proximity is unstable, the Approach 
Phase of the replenishment maneuver being of particular 
interest, where, excluding human intervention, collision is 
unavoidable. 

The mathematical modeling of two ships as one system 
With multi-inputs and multi-outputs has proven to be a 
powerful technique in the study of the dynamic behaviour of 
vehicles coupled by a common media. 

The station keeping control systen, initially based on a 
previous concept, was tested under di£rerent conditions, 
and, significant modifications were introduced, mainly in 
the gain parameters of the feedback loops, to satisfy the 
aims of safety as well aS precision in the whole maneuver. 

Фо avoid increasing the order of the general systen 
quations, nc forward speed control was included, being 
assased that the forward speed would bs controlled manually 
s lua Ly dones 

Even though the design of the hardware was not intended, 
it is eStimated that a Nodular System composed of a Doppler 
Sonar, a Gyroscopic Element and probabiy a Digital Logic 


ШОП бе suitable in its realization. 


Y 








Figure Li-1. 


rientation of the Space Axis 


and the Moving Axis (X, Y). 
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